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Abstract 

The  terahertz  (THz)  electromagnetic  frequency  band  is  a  potentially  useful  domain  for  remote 
sensing  in  military  and  Homeland  Security  applications.  Sensors  operating  in  this  band  gain 
some  of  the  material  penetrating  abilities  of  radio  waves,  some  of  the  image  resolving 
capabilities  of  infrared  (IR)  devices,  and  an  enhanced  chemical  discrimination  capacity  inherent 
to  the  THz  band  itself.  Both  passive  (thermal)  and  active  (maser)  THz  sensor  operating  modes 
have  been  contemplated,  with  the  greatest  emphasis  on  active  sensors.  Active  sensors  gain  a 
particular  advantage  due  to  typical  equivalent  source  brightness  temperatures  in  the 
neighborhood  of  1018  K.  At  these  intensities,  the  significance  of  natural  emission  effects  is 
negligible,  but  the  influence  of  water  vapor  absorption  can  be  significant  at  frequencies  above 
0.4  THz.  For  systems  performance  analysis  we  have  developed  simulation  software  to  assess  the 
limiting  effects  of  the  atmosphere  on  both  passive  and  active  terahertz  imagers  for  various 
instrumental  operating  frequencies,  bandwidths,  sensor  noise  levels,  integration  times,  and 
atmospheric  conditions.  We  will  show  some  of  the  early  results  of  our  modeling  and  simulation 
efforts,  which  illustrate  the  wide  range  of  potential  for  atmospheric  impacts  on  terahertz  sensors. 

1.  Introduction 

Recent  interest  in  the  terahertz  region  of  the  electromagnetic  spectrum  (arbitrarily  defined  as 
frequencies  in  the  range  0.4-10  THz)  has  prompted  development  of  imaging  systems  for  this 
band  in  both  the  astronomical  and  terrestrial  remote  sensing  communities.  For  the  latter  group, 
the  advantage  for  sensors  operating  in  the  THz  band  is  that  they  have  the  potential  to  combine 
the  penetrating  ability  of  radio  frequency  radiation  with  the  sensitivity  to  chemical  composition 
that  is  characteristic  of  infrared  (IR)  wavelengths.  Thus,  besides  potential  applications  such  as 
non-destructive  screening  of  packages  for  hazardous  materials,  THz  sensors  can  penetrate  layers 
of  clothing  to  detect  contraband.  These  operating  modes  can,  in  many  instances,  involve  indoor 
and  outdoor  operation  over  line-of-sight  (LOS)  distances  of  tens  or  hundreds  of  meters. 

Unfortunately,  THz  band  radiation  (sometimes  termed  T-rays)  is  strongly  absorbed  by 
atmospheric  water  vapor,  making  it  very  sensitive  to  ambient  absolute  humidity,  which  restricts 
the  range  of  practical  detection  for  passive  and  (to  a  lesser  degree)  active  THz  sensors. 
Absorption  and  scattering  of  THz  radiation  may  also  be  significant  for  heavy  rain,  snow,  and  fog, 
due  to  the  large  size  of  these  atmospheric  particulates  relative  to  THz  wavelengths  (which  are  on 
the  order  of  0.5  mm). 
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As  part  of  the  Defense  Advanced  Research  Projects  Agency  (DARPA)  funded  effort, 
coordinated  by  the  U.S.  Army  Night  Vision  and  Electronic  Sensors  Directorate  (NVESD),  the 
U.S.  Army  Research  Laboratory  (ARL)  Battlefield  Environment  Division  is  developing  system 
performance  and  atmospheric  effects  models  for  generic  THz  sensors.  In  this  paper,  we  shall 
discuss  that  portion  of  the  ARL  program  concerned  with  the  characterization  of  atmospheric 
effects  upon  prospective  THz  imaging  systems. 

2.  General  Simulation  Approach 

The  basis  of  our  image  propagation  approach  involves  simulation  of  a  pinhole-type  camera 
including  an  image  plane,  system  lens,  and  objective  field.  Within  this  simulated  space,  pixels  in 
the  image  plane  are  mapped  as  they  pass  through  the  object  space  until  they  reach  a  termination 
point  in  the  external  environment.  This  backward  ray  tracing  method  sums  incremental  path 
segment  contributions  for  each  ray  to  both  the  total  atmospheric  path  radiance  and  the 
transmitted  background  scene  radiance.  The  terminus  of  a  ray  trace  occurs  when  the  trace  hits  a 
target  or  background  surface  or  the  cumulative  optical  depth  exceeds  a  very  large  preset  value. 
Once  this  maximum  value  is  reached,  further  significant  incremental  contributions  to  the 
cumulative  path  radiance  are  unexpected.  This  cutoff  also  avoids  unnecessary  computations  for 
trace  directions  that  traverse  open  atmosphere  and  will  not  intersect  any  solid  target  or 
background  surface. 

The  model  sensor  (Figure  1)  in  the  initial  versions  of  our  simulation  was  a  simple  large  F-ratio 
camera  of  focal  length / and  with  an  N  x  N  square  array  of  uniformly-spaced  pixels  in  its  image 
plane.  The  physical  length  /  of  each  side  of  the  pixel  array  is  specified  by  the  user,  as  is  the  focal 
length  f  The  value  of  N  is  typically  set  to  either  256  or  5 12  pixels,  depending  on  run  speed 
versus  fidelity  tradeoffs.  Two  types  of  targets  are  modeled:  three-dimensional  (3D)  rectangular 
prism  targets  and  billboard  targets.  The  former  target  type  has  the  advantage  of  being  more 
physically  realistic  from  the  propagation  standpoint,  but  it  is  of  a  simplistic  structure.  The  latter 
class  of  targets  can  take  actual  THz  short  range  imagery  and  simulate  the  appearance  of  these 
targets  at  longer  ranges,  with  the  disadvantage  that  the  background  radiances  may  be  inaccurate, 
particularly  if  a  passive  billboard  target  is  moved  from  an  indoor  to  an  outdoor  situation. 


Figure  1.  Schematic  diagram  of  backward  ray  trace  geometry. 


The  atmospheric  environment  through  which  the  image  ray  tracing  occurs  is  modeled  as  a  3D 
gridded  space  in  which  the  amount  of  humidity  varies  from  position  to  position  in  discontinuous 
increments  via  a  system  of  volume  elements  (voxels).  The  length  of  each  voxel  is  set  by  the 
user.  Due  to  constraints  of  the  computing  environment,  the  voxel  pattern  is  modeled  using  a 
128  x  128  x  128  periodic  tiled  cubical  array. 

3.  Model  Development 

A  crucial  prerequisite  for  our  THz  image  propagation  model  is  an  adequate  statistical 
characterization  of  the  water  vapor  density  fluctuations  observed  in  the  atmosphere.  We 
approached  this  problem  through  a  series  of  field  measurements  using  a  fast  hygrometer  (model 
LI-7500,  developed  by  LI-COR®,  Inc.)  based  on  short-path  transmissometry  technology.  This 
instrument  was  used  to  sense  the  advected  water  vapor  concentration  field  at  a  20  Hz  sampling 
rate.  Sonic  anemometers  located  close  to  this  sensor  measured  the  wind  speed  and  provided  a 
means  to  convert  measurements  at  discrete  times  into  absolute  humidity  measurements  at 
positions  along  a  “wind  run”  axis.  Figure  2  shows  a  sample  of  the  data  collected  by  an  LI-7500 
and  a  nearby  sonic  anemometer.  Figure  3  shows  these  same  data  after  conversion  of  the 
temporal  axis  to  a  spatial  (wind  run)  axis,  and  rescaling  of  the  water  vapor  concentration  (in  units 
of  mmol  m~3)  data  into  an  absolute  humidity  measure  (in  units  of  g  m'3). 

The  field  measurement  data  were  then  analyzed  to  produce  a  one-dimensional  (ID)  power 
spectrum  of  the  absolute  humidity  fluctuations.  After  conversion  of  the  time-series  data  for 
absolute  humidity  to  the  wind  run  format,  a  fast  Fourier  transform  was  applied  to  compute  a  ID 
power  spectrum  for  each  one-hour  batch  of  data.  The  spectra  were  fairly  noisy,  so  a  5  point 
moving  average  was  used  to  smooth  them  so  that  their  structural  form  could  be  more  readily 
compared  to  the  known  spectral  behaviors  of  other  conservative  passive  additive  atmospheric 
properties  (e.g.,  Kaimal  et  al.,  1972).  As  expected,  the  spectra  showed  significant  departures 
from  power  law  behavior  at  very  low  and  very  high  spatial  frequencies,  where  source  and 
dissipation  phenomena  should  dominate,  respectively. 

The  intermediate  spatial  frequencies  (also  known  as  the  “inertial  subrange”)  did  show  conformity 
with  a  negative  5/3  power  law  spectrum  that  is  characteristic  of  a  ID  Kolmogorov  fluctuation 
spectrum.  The  3D  manifestation  of  this  spectrum  follows  a  negative  1 1/3  power  law  form. 

Figure  4  illustrates  the  appearance  of  the  raw  and  smoothed  spectrum  for  the  processed  data  of 
figures  2  and  3,  along  with  a  negative  5/3  power  law  fit  line. 


Figure  2.  Water  vapor  concentration  (upper  curve)  and  horizontal  wind  speed  (lower  curve) 
measured  for  a  1600  s  period  starting  at  1000  MST  27  March  2005  at  White  Sands 
Missile  Range  (WSMR),  NM. 
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Figure  3.  WSMR  water  vapor  concentration  data  for  1000  MST  27  March  2005  after 
conversion  to  absolute  humidity  and  wind  run  scaling. 
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Figure  4.  Spectral  power  density  derived  from  water  vapor  concentration  data  measured  over  an 
hour  period  starting  at  1000  MST  27  March  2005  at  WSMR. 

NOTE:  The  blue  curve  is  the  5  point  moving  average.  The  green  line  is  the  negative  5/3  power  law  curve. 


These  power  spectrum  results  are  employed  in  our  simulation  to  yield  a  stochastic  (fluctuating) 
absolute  humidity  field pr(t)  through  the  following  relation  (O’Brien  and  Tofsted,  2006): 


Pr  (r)  =  A) 


WII 


Nr(S)^\%)  exp(z' 2x  r  •%)  d\ 


(1) 


where  r  is  the  spatial  position  vector;  p0  is  the  spatial  average  of  absolute  humidity;  %  is  the 
spatial  frequency  vector;  Nr(£)  is  a  complex  Gaussian  random  variate  that  has  unit  variance  and  a 
zero  mean  in  both  its  real  and  complex  parts;  and  vPp(^)  is  currently  the  negative  1 1/3  power  law 
(Kolmogorov)  3D  spectrum.  The  second  term  on  the  right  hand  side  of  equation  (1)  is  the 
inverse  3D  Fourier  transform  of  the  function  Nr(Q  (vFj0(^))1/2.  The  square  root  arises  from  the  3D 
Fourier  transform  relationship  between  the  power  spectrum  and  the  covariance  of  the  absolute 
humidity  fluctuations.  When  a  Kolmogorov  spectrum  was  fitted  to  and  differenced  with  the 
observed  spectral  data  (the  zero  mean),  unit  variance  behavior  was  found  to  be  a  reasonable 
characterization  for  the  inertial  sub-range.  The  imaging  model  uses  discrete  formulations  of  the 
above  expressions  for  'P P(Q,  Nr(^),  and  pn(r),  and  a  discrete  FFT-based  inverse  transform 
(Tofsted,  2001). 


Another  critical  component  for  our  simulation  is  the  1993  version  of  the  Millimeter- wave 
Propagation  Model  (MPM93)  (Liebe  et  al.,  1992  and  1993).  This  model  enables  us  to  convert 
the  absolute  humidity  and  temperature  fields  into  a  stochastic  extinction  coefficient  field  at 
various  THz  frequencies,  which  are  used  to  directly  compute  ray  trace  radiances  and 
transmissions  in  our  imaging  model.  A  typical  MPM93  result  for  absorption  coefficients 
spanning  the  0.1  to  1  THz  frequency  band  is  shown  in  Figure  5.  Our  imaging  model 


encapsulates  MPM93  results  over  a  wide  range  of  frequencies,  atmospheric  temperatures,  and 
relative  humidities  in  the  form  of  a  look-up  table  (LUT)  that  is  interpolated  to  provide  THz 
extinction  coefficients  for  arbitrary  absolute  humidity,  temperature,  and  frequency. 


- 50%  Relative  Humidity  - Dry  Air 


Figure  5.  Comparison  of  MPM93  absorption  coefficients  for  moist  (50%  relative  humidity)  and  dry  (0% 
relative  humidity)  air  at  15  °C  and  1  atmosphere. 


Once  a  stochastic  absolute  humidity  field  has  been  computed  through  application  of  equation  (1) 
and  the  MPM93  model,  the  model  proceeds  to  ray  trace  from  each  pixel  position  in  the  image 
plane  back  to  object  points  in  the  external  environment.  For  passive  calculations,  where  no 
active  THz  sources  illuminate  the  target,  the  total  spectral  radiance  L(s )  for  an  LOS  path  of 
length  5  is  the  sum  of  the  attenuated  target  radiance  and  incremental  path  radiance  over  a  path  of 
N  segments  (Hoock  et  al.,  1992;  O’Brien  and  Tofsted,  1998  and  2006): 


ZAT.As,  [l-expf-KA.Ik,, 


j= 0 


(2) 


where  Lj{s)  is  the  unattenuated  target  radiance;  Aj  and  A,v,  are  the  respective  extinction 
coefficient  and  length  of  line  segment  i;  and  LSi  is  the  limiting  path  radiance.  The  limiting  path 

radiance  is  the  radiance  that  would  result  as  the  optical  depth  of  the  medium  tends  to  large 
values.  If  scattering  is  neglected  (which  is  a  good  approximation  in  the  THz  band  in  the  absence 
of  heavy  rain,  snow,  or  fog),  the  limiting  path  radiance  is  essentially  the  blackbody  radiance  at 
ambient  atmospheric  temperature.  Our  imaging  model  currently  follows  this  assumption,  since  it 
does  not  yet  incorporate  THz  hydrometeor  scattering  algorithms. 


Three  variants  of  our  imaging  model  exist  to  handle  different  sensing  modes:  1)  spectral  (single 
frequency)  passive,  2)  broadband  passive,  and  3)  active.  The  spectral  and  broadband  modes  use 
equation  (2)  to  calculate  pixel  radiances.  The  broadband  mode  sums  over  frequency  bins 
weighted  by  sensor  response  and  includes  sensor  sensitivity  via  a  noise  equivalent  power  (NEP) 
input.  Thus,  even  with  contrast  stretching,  sensor  noise  will  limit  the  target  contrast  once  the 
attenuated  target  signal  falls  below  a  threshold  level.  Examples  of  this  phenomenon  will  be  seen 
in  section  4.  The  active  mode  assumes  that  a  narrow  band  source  illuminates  the  target  and  that 
only  the  first  term  of  equation  (2),  multiplied  by  an  outgoing  transmission  factor,  is  a  significant 
contributor  to  the  apparent  pixel  radiance.  This  assumption  is  motivated  by  the  fact  that  even 
modest  active  THz  sources  have  enormous  effective  blackbody  temperatures  in  the  THz  region. 
For  the  active  imaging  examples  shown  below,  we  used  the  nominal  source  characteristics  for 
the  active  source  used  by  the  Ohio  State  University  (OSU)  group  in  laboratory  imaging,  which 
had  an  effective  source  brightness  temperature  of  1018  K. 

4.  Results 

We  first  exercised  the  passive  spectral  radiance  version  of  our  model  for  3D  targets  with  uniform 
or  simple  periodic  surface  emissivities.  Figure  6  illustrates  how  the  appearance  of  a  rectangular 
box  target  might  vary  as  its  surface  temperature  under  fixed  atmospheric  conditions.  We  have 
used  MPM93  results  to  establish  a  mean  value  for  the  extinction  at  a  study  frequency  of  650 
GHz.  As  can  be  seen  from  figure  5,  650  GHz  is  a  frequency  where  water  vapor  extinction  is 
near  a  local  minimum  or  “window.”  The  results  shown  in  figure  6  use  the  Kolmogorov  spectrum 
for  fluctuations  about  this  mean,  as  implemented  from  equation  (1).  The  ground  plane  is  at  a 
temperature  of  273  K,  the  atmospheric  mean  temperature  is  275  K,  and  both  the  box  and  ground 
plane  are  assigned  reflectivities  of  0.2.  A  vertical  bar  pattern  is  imposed  on  one  of  the  target  box 
faces,  with  bar  temperatures  alternating  between  the  nominal  target  value  and  that  of  the  ground 
plane. 

The  target  temperature  is  depressed  10  K  below  the  ground  plane  temperature  in  figure  6(a)  and 
displays  strong  negative  contrast  in  the  contrast-stretched  image.  The  same  holds  true  for  figure 
6(b),  where  the  target  temperature  is  depressed  by  5  K,  although  some  lightening  is  evident  on 
the  target  sides  and  bar  pattern.  The  target  is  at  the  same  temperature  as  the  ground  plane  in 
figure  6(c),  as  evidenced  by  the  disappearance  of  the  ground  temperature/target  temperature  bar 
pattern  and  the  indistinctness  of  the  edges  of  the  target  that  rest  on  the  ground  plane.  However, 
strong  negative  contrast  is  visible  on  all  of  the  elevated  target  edges  due  to  absorption  and 
emission  from  the  background  atmosphere.  Figures  6(d)  and  6(e)  represent  a  1  K  and  2  K 
elevation  of  target  temperature  above  that  of  the  ground  plane.  As  one  might  expect,  some 
vertical  target  edge  regions  begin  to  show  positive  contrast,  but  the  upper  target  edges  still 
display  appreciable  negative  contrast.  Only  with  the  target  temperature  elevated  5  K  over  the 
ground  plane  temperature  and  3  K  above  the  atmospheric  temperature  do  we  finally  see  near  zero 
contrast  on  the  upper  edges,  as  seen  in  figure  6(f). 


Figure  6.  Variation  of  650  GHz  target  contrast  with  target  temperature  in  a  Kolmogorov  spectrum  emission/absorption 
atmosphere.  Target  temperatures  were  (a)  263  K,  (b)  268  K,  (c)  273  K,  (d)  274  K,  (e)  275  K,  and  (f)  278  K. 
Ground  and  air  temperatures  were  273  K  and  275  K  respectively. 


Next,  we  considered  a  more  realistic  sensor  that  incorporated  operational  properties  of  passive 
scanned  bolometric  sensors  similar  to  those  designed  by  the  TIFT  team  at  OSU.  We  assumed  a 
nominal  passband  of  100  GHz  centered  at  a  frequency  of  650  GHz  and  integrated  in  25  GHz 
steps.  The  0.3  m  aperture  sensor  had  a  frame  rate  of  100  Hz  and  a  detector  NEP  of  3  x  10"12  W 
Hz  1/2,  yielding  a  noise  power  of  3  x  10'11  W.  The  random  additive  pixel  noise  in  the  detector 
was  then  computed  as  the  noise  power  multiplied  by  a  Gaussian  random  number  variate  of  zero 
mean  and  unity  standard  deviation.  This  noise  component  was  applied  to  ray  traces  over  a  60  x 
60  pixel  image  plane  for  a  billboard  target  that  is  a  THz  image  of  a  subject  with  a  mock-up  pistol 
concealed  beneath  a  layer  of  clothing.  The  results  of  the  simulation  are  displayed  in  figure  7  for 
two  different  propagation  ranges  (40  m  and  80  m)  and  four  levels  of  average  absolute  humidity 
(0.7,  7.6,  17.9,  and  35.7  g  m‘3).  The  results  clearly  demonstrate  some  of  the  shortcomings  of  the 
passive  imaging  mode.  At  80  m,  all  of  the  images  have  severe  noise  degradation  of  the  signal, 
even  for  the  relatively  cold  and  dry  scenario  on  the  left-hand  column  of  the  figure.  The  40  m 
cases  show  reasonable  performance  up  to  absolute  humidities  in  the  neighborhood  of  7.6  g  m"3, 
but  noise  effects  are  still  quite  apparent  even  at  low  humidities. 


T  =  -10  °C 
Rel.  Hum.  =  30% 
Abs.  Hum.  -  0.7  g/m3 


T  =  20  °C 
Rel.  Hum.  =  44% 
Abs.  Hum.  -  7.6  g/m3 


T  =  43  °C 
Rel.  Hum.  =  30% 
Abs.  Hum.  -  17.9  g/m3 


T  =  35  °C 
Rel.  Hum.  =  90% 
Abs.  Hum.  -  35.7  g/m3 


Figure  7.  Broadband  passive  imaging  examples  centered  at  650  GHz.  Images  in  the  upper  row  are  at  a  range  of 
40  m,  those  in  the  bottom  row  are  at  80  m. 


The  OSU  group  also  has  developed  short  range  active  THz  systems.  Figures  8  and  9  show 
results  for  a  simulated  propagated  active  mode  image  at  650  GHz,  for  sensor  dynamic  ranges  R 
(maximum  detectable  radiance  divided  by  minimum  detectable  radiance)  of  80  and  50  dB,  with  a 
threshold  minimum  detectable  power  Lo=  10'3  x  1  0’ft/|°  W.  The  source  power  was  assumed  to 
be  1  mW  for  the  examples  shown  here,  which  use  the  same  set  of  path  lengths  and  atmospheric 
conditions  as  the  passive  cases  above. 


80  dB  Dynamic  Range 


T  =  -1 0  °C  T  =  20  °C  T  =  43  °C 

Rel.  Hum.  =  30%  Rel.  Hum.  =  44%  Rel.  Hum.  =  30% 

Abs.  Hum.  =  0.7  g/m3  Abs.  Hum.  =  7.6  g/m3  Abs.  Hum.  =  17.9  g/m3 


T  =  35  °C 
Rel.  Hum.  =  90% 
Abs.  Hum.  =  35.7  g/m3 


Figure  8.  Active  imaging  example  for  650  GHz  source.  The  dynamic  range  of  sensor  was  set  to  80  dB.  The  top 
row  is  at  a  range  of  40  m,  the  bottom  row  is  at  80  m. 


50  dB  Dynamic  Range 


T  =  -10  °C  T  =  20  °C 

Rel.  Hum.  =  30%  Rel-  Hum.  =  44% 

Abs.  Hum.  =  0.7  g/m3  Abs.  Hum.  =  7.6  g/m3 


T  =  43  °C 
Rel.  Hum.  =  30% 
Abs.  Hum.  =  17.9  g/m3 


T  =  35  °C 
Rel.  Hum.  =  90% 
Abs.  Hum.  =  35.7  g/m3 


Figure  9.  Active  imaging  example  for  650  GHz  source.  The  dynamic  range  of  sensor  was  set  to  50  dB.  The  top 
row  is  at  a  range  of  40  m,  the  bottom  row  is  at  80  m. 


The  active  imaging  mode  shows  much  better  performance  than  does  the  passive  mode  through 
equivalent  atmospheric  conditions.  With  an  80  dB  sensor  dynamic  range,  target  features  are 
significantly  degraded  only  at  the  longest  range  and  highest  absolute  humidity.  The  imaging 
performance  degrades  at  considerably  lower  absolute  humidities  for  the  50  dB  case,  however. 


5.  Conclusions  and  Future  Work 


The  THz  imaging  model  described  in  this  paper  shows  promise  as  a  means  to  investigate 
atmospheric  impacts  on  both  real  and  potential  imaging  systems.  However,  there  are  several 
areas  where  the  package  must  be  significantly  improved  before  this  capability  becomes  a  reality. 
First,  the  characterization  of  the  geometrical,  optical,  and  thermal  properties  of  the  3D 
environment  must  be  made  more  complete.  Existing  target  models  employed  by  target 
acquisition  packages  need  to  be  investigated  as  resources  that  might  satisfy  part  of  this 
requirement.  Another  significant  part  of  the  THz  environmental  description  task  might  be 
addressed  by  studies  of  natural  and  manmade  background  surface  emissivity  and  thermal 
characteristics  in  the  millimeter  wave  literature. 

We  will  also  need  to  improve  the  sensor  model  to  incorporate  the  more  detailed  real  world 
properties  of  optical  and  electronic  systems  that  we  know  will  significantly  impact  performance. 
In  particular,  we  will  need  to  apply  a  modulation  transfer  function  (MTF)  methodology  to  the 
image  formation  process.  The  engineering  models  that  characterize  these  imperfections  of  the 
sensor  system  are  well  known,  but  as  yet  we  do  not  have  well  defined  real-world  THz  system 
designs  to  which  they  may  be  applied.  Differences  in  image  appearance  will  also  depend  on  the 
scanning  technique  used.  Instantaneous  scanning  of  all  pixels  produces  one  result.  Single-pixel 
scanning  allows  the  humidity  field  to  evolve  during  the  image  acquisition  process,  etc. 

If  hydrometeors  are  present  in  the  atmosphere  for  operational  scenarios,  it  will  also  be  very 
important  to  add  a  hydrometeor  scattering-absorption  algorithm  to  the  active  source  model. 
Liebe’s  MPM93  (or  a  similar)  model  provides  the  absorption  component  for  snow  and  rain 
attenuation  of  the  source  beam,  but  it  does  not  directly  supply  any  estimate  for  the  scattered 
radiance.  This  effect  could  significantly  modify  the  path  radiance  signal.  In  the  absence  of  such 
resonance-scale  objects  as  snow  crystals  or  rain  drops,  the  atmospheric  water  vapor  will  simply 
attenuate  the  illumination  beam.  The  terahertz  reflectance  properties  for  target  and  background 
surfaces  will  also  be  critical  scenario  elements. 
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Introduction  - 1 


•  Key  Background  Points: 

-  Recent  advances  in  terahertz  imaging  sensor 
technology  are  leading  to  higher  frame  rates 
and  pixel  density 

-  Some  terahertz  sensor  applications  will 
involve  outdoor  operation  in  various  climatic 
environments 

-  Atmospheric  absorption  is  strong  in  the  low 
terahertz  band  (0.5  -  1  GHz),  and  is  due 
primarily  to  water  vapor 


Introduction  -  II 


•  Issues  for  Outdoor  Terahertz  Operation: 

-  Statistics  for  spatial  variation  of  water  vapor  mass 
density  (or  absolute  humidity)  about  the  mean  are 
poorly  known,  particularly  as  functions  of  terrain  and 
prevailing  weather  conditions 

-  Functional  form  for  the  spectral  power  density  of 
absolute  humidity  is  not  well  known 

-  Target  contrast  as  a  function  of  range  for  terahertz 
active  and  passive  imagers  is  not  well  established 


Introduction  -  III 


•  Simulation  Approach: 

-  Collect  preliminary  water  vapor  fluctuation  field  data 
and  create  analysis  methods  and  software 

-  Analyze  absolute  humidity  data  to  establish 
reasonable  functional  shape  for  spectrum  and 
amplitude  for  fluctuations 

-  Use  readily  available  models  to  connect  simulated 
stochastic  absolute  humidity  and  air  temperature 
fields  with  absorption  and  radiance  fields  at  terahertz 
frequencies 

-  Implement  a  ray-trace  imaging  model  for  billboard 
and  three-dimensional  targets  in  the  simulated 
stochastic  atmosphere 


MPM93  Absorption  Coefficient 
Results:  TAtm  =  25C,  35C 
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MPM93  Sample  Results:  25  °C 
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WSMR  Water  Vapor 
Measurements 


Measurement  Periods:  23-28  Mar  05,  late  Aug  05 
Measurement  Sampling 

-  20  Hz  rate 

-  1  hour  data  batches 

-  Continuous  over  day 
Measurement  Apparatus 

-  Single  Li-Cor  Model  LI-7500  absolute  humidity  sensor 
system 

-  Four  sonic  anemometer  probes  and  support 
equipment 


Example  from  August  05 
Las  Cruces  Measurements 
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WSMR  Water  Vapor  Measurement 
Results:  23  Mar  05  1600  MST 


WSMR  Water  Vapor  Measurement 
Results:  26  Mar  05  1 000  MST  W 
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Absolute  Humidity  Variations 
on  Different  Time  Scales 

1  Second  Averages  of  20  Hz  Data 
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WSMR  Water  Vapor  Measurement 
Results:  26  Mar  05  1000  MST 
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Higher  Absolute  Humidity  Example 


Temperature  Spectrum  Example 


1400  -  1600,  25  Aug  2005,  6’0”  Hgt 
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Simulation  of  Atmospheric  Fluctuations 


•  The  one-dimensional  distribution  of  absolute  humidity  p(x)  may  be 
expressed  as  the  sum  of  an  average  />0(x)  and  a  randomly  fluctuating 
component  />i(x): 


p(x)  =  p0(x)  +  p1(x) 

•  The  randomly  fluctuating  component  has  an  expectation  value  of 
zero: 

(Pi(x))  =  0 

•  Under  the  assumption  that  the  spatial  expectations  and  spatial 
averages  are  equivalent  (homogeneity),  the  correlation  of  p  values  at 
two  different  locations  Xi  and  x2  may  be  expressed  using  the 
covariance  function 

rp(Xi,x2)  =  <  [p(x1)-<p(x1)>]  [p(x2)-<p(x2)>]  >, 

or 

Tp(x)  =  <  p1(x1)p1(x1+x)> 


Simulation  of  Atmospheric  Fluctuations 


•  A  three-dimensional  covariance  value  for  the  random  fluctuation  pR 
may  then  be  defined  by 


00 

rKR(r)  =  j  J  I  Nr(9  Vk(A)  e12””^ 


•  This  result  is  then  be  used  to  define  an  expression  for  the  random 

fluctuation  field  about  a  mean  value: 

00 

PR(r)  =  Po+j  J  J  Nr(©T“(9  el2””? 


where  the  three-dimensional  Kolmogorov  inertial  sub-range  spectrum 
with  integration  constant  C  is  given  by 


x¥(Q  =  C^n/3 


with  the  spectrum  and  covariance  related  by  the  transform  pair 
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Synthetic  Image  Example 

•  650  GHz,  albedo  =  0.2  for  target  and  ground 

•  2  m  x  1 .5  m  x  0.6  m  box  target 

•  -X  face  of  target  has  alternating  Tgnd/Ttgt 
vertical  stripes 

•  Sensor  at  (x,y,z)  =  (0,  0,  15  m) 

•  Target  at  (50  m,  40  m,  0.3  m) 

•  Sensor  focal  length  150  mm,  10  mm  square 
focal  plane,  512  x  512  pixels 

•  Tatm,  Tgnd  fixed  at  288  K  (15  °C),  285  K 

•  100%  RH,  Kolmogorov  spectrum,  20%  a 


284  K  285  K  286  K 


Planned  Terahertz  Atmospheric  \  f 
Imaging  Model  Upgrades 

•  Active  imaging 

-  Bistatic  and  monostatic  options 

•  Hydrometeor  treatment 

-  Rain  and  snow  absorption  and  scattering 

•  Simplified  user  input 

-  Relative  humidity  instead  of  absolute  humidity 

-  Option  for  user  input  of  temperature  and  humidity  variances 

•  More  choices  for  sensor  performance  characteristics 

-  Focal  plane  array  or  raster  scan  image  formation  methods 

-  More  realistic  system  noise  and  angular  resolution  properties 

•  Obtain  better  terahertz  band  background  imagery  data 
as  it  becomes  available 


Planned  Field  Support  Work 

•  Continue  analysis  of  WSMR/Las  Cruces  Li-Cor 
and  sonic  anemometer  data 

•  Collect  and  analyze  WSMR  data  sets  using  two 
Li-Cor  sensors  at  variable  spacings  to  get 
absolute  humidity  covariance 

•  Collect  and  analyze  concurrent  Li-Cor  and  sonic 
anemometer  data  sets  to  determine  cross 
variance  between  temperature  and  absolute 
humidity 
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WSMR  Water  Vapor  Measurement 
Results:  27  Mar  05  1000  MST 
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WSMR  Water  Vapor  Measurement 
Results:  27  Mar  05  1000  MST 


MPM93  Sample  Results:  T  =  25°C, 
Frequency  =  600  GHz 
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